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The magnetic properties of the two-channel periodic Anderson model for uranium ions, comprised 
of a quadrupolar and a magnetic doublet are investigated through the crossover from the mixed- 
valent to the stable moment regime using dynamical mean field theory. In the mixed-valent regime 
ferromagnetism is found for low carrier concentration on a hyper-cubic lattice. The Kondo regime 
is governed by band magnetism with small effective moments and an ordering vector q close to 
the perfect nesting vector. In the stable moment regime nearest neighbour anti-ferromagnetism 
dominates for less than half band filling and a spin density wave transition for larger than half 
filling. T m is governed by the renormalized RKKY energy scale Me// <X„> 2 J 2 po(aO- 
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Introduction. Heavy Fermion (HF) materials [EJ have 
drawn much attention since the discovery of supercon- 
ductivity in CeCu2Si2 ^| which is characterized by an 
anisotropic order parameter with symmetry yet to be de- 
termined. Over the years, it became clear that the heavy 
Fermi-liquid is not the ground state: almost all HF ma- 
terials show magnetic or superconducting phase transi- 
tions, which cither compete with each other in Ce based 
compounds or coexist in Uranium based materials. Addi- 
tionally, doubts were casted on whether the Fermi-liquid 
picture is even applicable for the normal state properties 
in some U-based || or doped compounds Q. 

In this paper, new results on magnetism in the two- 
channel periodic Anderson model (TCPAM) || are pre- 
sented. In the TCPAM, the restricted local Hilbert space, 
consisting of a quadrupolar and a magnetic ground state 
doublet for two different valence configurations, models 
an Uranium ion in a cubic environment, which can either 
fluctuate between 5/ 2 /5/ 3 or 5/ 1 /5/ 2 ground state con- 
figurations. On symmetry grounds, the coupling of the 
ion to the itinerant electrons via hybridization can be 
characterized by an orbital and a spin label generating 
an effective two band model |J. 

It has been shown that the paramagnetic phase ex- 
hibits non-Fermi-liquid properties, as manifested in a 
large residual resistivity [0 . They are caused by the over- 
compensated effective impurity on each Uranium site act- 
ing as incoherent scatterer. The predicted absence of a 
Drude peak in the optical conductivity has been found 
experimentally in UBei3 [||. On a lattice, however, the 
residual entropy @] is expected to be quenched by a phase 
transition. Candidates are magnetic, channel ordering 
or superconducting transitions, which all have to bind 
the residual magnetic or quadrupolar spin. In the peri- 
odic two-channel Kondo model, anti-ferromagnetic and 
charge ordering phase transitions have been reported in 
the strong coupling limit of the model |J|po(0) > 0.5 [p|, 



Po(e) being the non-interacting density of states. Addi- 
tionally, the authors interpret a sign change in the pair 
susceptibility for a spin/channel singlet sector as a su- 
perconducting transition. The absence of a divergence is 
viewed as indicator of a first order transition. 

In this paper, I focus on possible spin transitions. I 
have calculated the spin phase transition temperature 
T m as a function of band filling and energy difference 
f-f = E a — E a between the quadrupolar and the magnetic 
doublet using the dynamical mean field theory |^] . Vary- 
ing e/ continuously allows me to investigate three differ- 
ent regimes: |e/| smaller than the charge fluctuation en- 
ergy A defines the intermediate valence regime, which has 
recently been discussed experimentally [|f0],[ll| and theo- 
retically ffj,[]l| as a possible scenario for Ui-^Th^Beia. 
The local moments are not well defined, and the carriers 
rather mobile. Magnetic ordering would have to involve 
a gain in kinetic energy since the effective exchange cou- 
plings are small p| ]. Ferromagnetism is found at low 
carrier concentration. For |e//A| « 2, the Kondo effect 
dominates. Band magnetism with very small effective 
magnetic moments is expected since all available elec- 
trons are involved in the screening process. For large 
|e//A|, magnetism is governed by stable local moments 
whose interactions are mediated by the conduction elec- 
trons (RKKY interaction). 

Theory. The two channel Anderson lattice Hamilto- 
nian under investigation reads 



E^R^L + ^} • (i) 



a<i,j> 



X are the Hubbard operators, d being the spatial dimen- 
sion, and t* the reduced hopping matrix element of the 
conduction electrons between nearest neighbours. The c- 
electrons are created by c\ aa at the lattice site i, labeled 
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by a spin a and a channel index a = ± and couple via 
the hybridization matrix element V to the ionic many- 
body states on each lattice site. We will assume that the 
energies of the magnetic doublets split linearly under an 
applied small external field, and van Vleck contributions 
are neglected. 

In the limit of infinite spatial dimensions the two 
particle- hole irreducible vertex Tff becomes local [ p"5| . 
Using the dynamical mean field theory (DMFT) in fi- 
nite dimensions, the Bethe-Salpeter equation for the full 
two particle-hole propagator in f-electrons reads 

a // (-«)=%° // +4^ // a // (-«) . (2) 

where x°ff is the free particle-hole //-propagator 

X°f f (iuj n ,iuj m ;w n ,q) = - 5 -^Yjk F k{iuJ n + iv n )Fk +1 (iu> n ) 
= J de 1 de 2 p°(e 1 ,e 2 ,q)F(ei,iu> n + M/„)F(e 2 , iw„) 

(3) 

and Fk(z) = F(e(k),z) is the lattice /-Green's func- 
tion. The matrices are indexed by Matsubara frequen- 
cies, Hff(iiu n ,iuj m ;w n ) = U 



ff 



i, an d p 3 ) (ei,e 2 ,q) 

is the two particle density of states given by the under- 
lying lattice Jl(J. Since Eqn. (0) also holds for the local 



two-particle propagator H^, we obtain the irreducible 
vertex Tff and substitute it into (^|) to derive 
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which is valid for each Boson frequency iv n separately 
pl| , p2[ . In order to study possible phase transitions, I re- 
stricted my calculations to the fully symmetric phase. I 
used a Gaussian non-interacting density of states po(e) = 
_i_ e -( e /* ) 2 f or a simple cubic lattice in large dimen- 
sions, and focused on the static //-spin susceptibility 



*//(£) = 7^2 



n // (iw n) icj m ;0, g )e 5 ^» +i — ) . (5) 



The diverges of Xffio) originating from zeros of the de- 
nominator of (^J) determine the phase transitions, resem- 
bling a Stoner like expression for the lattice susceptibil- 
ity. I would like to mention that the convergence factors 
in Eqn.Q are important for evaluating the susceptibil- 
ity analytically and numerically. The conduction elec- 
tron susceptibility Xcc and the mixed terms Xcf and Xfc 
also contribute to the total susceptibility. However, since 
the unperturbed conduction electrons form a Fermi-gas, 
any irreducible interaction is generated by the /-electrons 
[P0[| . Therefore, divergences occur collectively for all con- 
tributions when Xff(o) diverges due to c — / mixing. 



In large spatial dimensions, /^(ei, £2, q) obtains its q de- 
pendence only through lattice dependent structure factor 
n(q) = cos(qi)/D in the case of a hyper cubic lattice 
of £)-dimensions |l6|] . 

In context of Heavy- Fermion systems Q, the terms 
anti-ferromagnetism (AF) and spin density wave (SDW) 
are often used interchangeably, since it is hard to 
distinguish experimentally true incommensurable band 
magnetism from complicated yet commensurable mag- 
netic structures. Here, I will call a transition anti- 
ferromagnetic if the divergence occurs at n(q) = — 1, 
ferromagnetic, if n(q) = 1, while all other cases will be 
labeled as SDW. To be specific, I choose the wave vec- 
tor q = qo(l, 1, • ••) along the nesting direction instead of 
n(q) in all figures. I define the characteristic temperature 
of the lattice T* as the temperature at which the initial 
moment of the lower lying doublet is screened down to 
40% |0|. This energy scale is the lattice analog to the 
single ion Kondo temperature. Taking into account the 
lack of universality due to the ongoing renormalization 
of the conduction band for T < T* , it is not surprising 
if T* does not coincide with the energy scale of the elec- 
trical response Tfl, as demonstrated in Ref. (l?]]. Hints 
are also seen in calculations for low carrier concentrations 

Results. The Anderson width A = V r2 7rpo(0) is used 
throughout the paper as the energy unit. All calcula- 
tions have been performed for a fixed band width t* = 10 
which yields the hybridization strength V 2 — At* / If 
one defines the effective dimensionless coupling constant 
g = p (0)|J| = p Q (0)V 2 /\e f \ = A/( 7 r|e / |), and requires, 
that A/e/| < 1 as a criterion for excluding the inter- 
mediate valence regime, the TCPAM can only describe 
the weak coupling regime (g < l/V) of the two-channel 
Kondo lattice model. The chemical potential ijl has been 
varied from —5 to 5 to obtain different conduction elec- 
tron occupancies n c : at n c = 4 the c-bands are com- 
pletely filled. 

The occupancy of the magnetic doublet decreases 
monotonically with increasing c-band filling. Even 
though the effective DOS is altered symmetrically around 
half filling, the different /-valent states \a> (\a>) involve 
admixture of virtual hole (particle) band states which are 
increased (decreased) with decreasing chemical potential. 
Therefore, T* decreases more rapidly with conduction 
band depletion than with conduction band filling. More- 
over, assuming T* cx exp— l/(Jp), the change of T* be- 
comes more significant for decreasing J. This is exactly 
what is found numerically of the occupancy dependence 
of T*(n c ) for different e/. 

Generally, magnetic phase transitions in the periodic 
Anderson model are expected to be AF for the stable mo- 
ment or SDW for the Kondo regime W stemming from 
inter and intra band excitations in a hybridized band 
picture. This picture is not truly applicable in the two 
channel model due to lack of quasi-particles at the chem- 
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FIG. 1. Mixed- valent regime: T m /T* versus 50 for differ- 
ent conduction electron fillings. A ferromagnetic transition is 
found from 1/4 to 1/2 conduction band filling for e/ = — 1A. 
The inset shows T m /T* for e/ = 0. 

ical potential. 

In the mixed valent (IV) regime ( e//A = —1) the 
largest T m is found for a ferromagnetic transition at less 
than half band filling (Fig. 0). With the exception of 
n c = 1, T m /T* < 1. Independently, I also verified the 
existence of a ferromagnetic phase by calculating the spin 
polarization in an external magnetic field: its magnitude 
was decreased from a finite value of the order of 10 _4 T* 
to zero field. A zero field magnetization was found only 
below T m . Since the used XNCA method works reliably 
only in the multi-channel case , accurate determination 
of a magnetization curve below T m was not possible. 

In the single channel PAM, no magnetism is expected 
in the IV regime due weakness of the effective exchange 
couplings. In the two channel PAM, however, incoherent 
scattering of the conduction electrons remains present in 
the paramagnetic phase down to low temperatures. En- 
ergy can be gained by reducing the scattering rate in a 
spin-polarized state. This gain grows with increasing oc- 
cupancy of the magnetic doublet. Therefore, the transi- 
tion is found at low c-carrier concentration, where charge 
fluctuations mixes increasing number of c-hole states to 
local magnetic ground state, stabilizing the local mo- 
ment. Because charge fluctuations are ignored in the 
Kondo lattice model, it has not been observed before 
P). The FM is clearly a strong coupling effect, since 
it is caused by high energy charge fluctuations JlSp . It 
leads to a higher mobility of the c-electrons measurable 
through an enhancement of DOS in the majority spin 
above the chemical potential. In the FM phase, the in- 
creasing magnetization quenches the spin channel and 
the model crosses over to a single channel PAM || . Be- 
low T m , Fermi-liquid properties are regained. 

The picture changes significantly in the crossover 
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FIG. 2. Kondo regime: T m /T* versus qo for different 
conduction electron fillings, e/ = —2A. The inset shows 
Tm/T* near quarter filling. 

regime, e//A = —2 (Fig. ||): although the occupancy 
changes moderately from 0.82 to 0.74 with increasing 
band filling, the maximum of T m is now found at a fi- 
nite q vector qq « .85tt indicating a SDW transitions 
for band fillings larger than 0.5. Around half filling, the 
Kondo effect dominates. The effective local moment is 
already quenched to 7% of its high temperature value 
at the SDW transition temperature T m . The transition 
occurs well below T* while with increasing effective mo- 
ment, the transition is shifted towards higher tempera- 
tures. Near quarter filling, the local moments are only 
weakly screened: the effective moment ^// = 0-6^ ree , 
the ratio T m /T* is very large (inset of Fig. ||) and an 
AF transition is observed similar to the stable moment 
situation discussed next. 
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FIG. 3. Stable moment regime: T m /A versus qo for differ- 
ent conduction electron fillings, e/ = — 3.5A. 

Lowering e/ further, we arrive at the stable local mo- 
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ment regime (Fig. g). For n c < 2, an AF transition is 
found. Note that T m is not rescaled with T* in this case, 
indicating that the screening of the moments is very low 
(about 25%). Increasing the band filling, the admixture 
of \a> states into the ground state increases, which de- 
creases the local moments: for n c > 2 SDW transitions 
at large q vectors are found. For positive e/, no magnetic 
phase transition is found. 

Discussion. I have established, for the first time, a 
ferromagnetic solution of the two-channel periodic An- 
derson model in the intermediate valence regime for low 
carrier concentration. Energy and mobility is gained by 
reducing the incoherent scattering of the conduction elec- 
trons in a spin-polarized state. The loss of energy due to 
spin polarization of the conduction band is low for small 
band fillings. The transition is driven by high energy 
charge fluctuations indicating the strong coupling nature 
of the FM transition. The presences of FM transitions 
and strong ferromagnetic fluctuations suggest the pos- 
sibility of triplet superconductivity. Preliminary stud- 
ies indicate a divergence of the pair susceptibility in the 
spin/channel triplet sector, implying an odd-frequency 
gap function. 

In the Kondo regime, I find SDW transitions, with a 
large wave vector q, for band fillings larger than 1/2. 
Additionally, nearest neighbour AF and tendency to- 
wards ferromagnetism is observed for low carrier concen- 
trations. In the stable moment regime AF transitions are 
found at less than half filling and SDW transition above 
half filling. However, I have excluded a detailed treat- 
ment of band-structure effects in finite dimensions within 
DMFT, which can be included through Equ. (||). The 
bare two-particle density of states then becomes truly 
q-dependent. I hope to clarify its influence on the spin 
susceptibility in three dimensions in the near future. 

In all the cases I investigated, I found T m cx Erkky = 
nlfj <X a > 2 J 2 po(fi). where <X a > is the occupation 
number of the magnetic ground state doublet of the each 
ion. The ratio Erkky /T m decreases monotonically from 
around 3. to 1.5 with decreasing effective local moment 
for 6f = —2. The magnetic transition is always governed 
by an effective RKKY energy Erkky, which takes into 
account the screening of the local moments. Erkky can 
be one order of magnitude smaller than the bare RKKY 
energy. 

Experimentally, the change of parameters in real mate- 
rials is much more difficult: usually either the size of the 
unit cell is changed through isostructural substitution or 
through external pressure. Essentially, the overlap be- 
tween nearest neighbours is varied which yields a change 
in the hybridization V and the bandwidth t*. It is be- 
lieved, that the core energy difference et stays unaltered 
p3[ . Under the assumption, that only hybridization and 
bandwidth are altered by pressure, application of pos- 
itive pressure will be equivalent to a reduction of |e/|. 
The results presented show clearly the same tendency as 



observed in experiments |17|: the application of pressure 
leads to a suppression of magnetic phase transitions. Lo- 
cal moments are destroyed through hybridization induced 
de-localization of the f-clcctrons. 
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